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We study the decoherence properties of an entangled bipartite qubit system, represented by two 
two-level atoms that are individually coupled to non-Markovian reservoirs. This coupling ensures 
that the dynamical equations of the atoms can be treated independently. The non-Markovian 
reservoirs are described by a model which leads to an exact non-Markovian master equation of the 
Nakajima-Zwanzig form [J. Salo, S M. Barnett, and S. Stenholm, Opt. Commun. 259, 772 (2006)]. 
We consider the evolution of the entanglement of a two-atom state that is initially completely 
entangled, quantified by its concurrence. Collapses and revivals in the concurrence, induced by 
the memory effects of the reservoir, govern the dynamics of the entangled quantum system. These 
collapses and revivals in the concurrence are a strong manifestation of the non-Markovian reservoir. 

PACS numbers: 03.65.Yz, 03.65.Ud, 03.67.-a, 03.67.Mn, 42.50.-p 



I. INTRODUCTION 

The dynamics of open quantum systems is usually con- 
sidered under the assumption that the physical systems 
of interest are sufficiently isolated from their environ- 
ment. This allows one to utilize certain approximations 
such as the weak coupling and the Markovian approx- 
imation [l|. The first approximation ensures that the 
interaction between the system and the environment is 
sufficiently weak so that the system is quasi closed. The 
latter one assumes that the characteristic times of the 
reservoir are much shorter than that of the system. This 
allows one to ignore memory effects from the environment 
and one arrives at time-local equations of motion for the 
reduced state of the system which are ususally expressed 
in terms of Gorini-Sudarshan-Kossakowski-Lindblad op- 
erators |3, Q- These operators generate a completely 
positive time evolution which conserve the physical prop- 
erties of all quantum states 4'] . 

Recent studies have shown the limits of the Marko- 
vian description of quantum computation and quantum 
error correction 0, H, Qi S S Hy; Hll ■ In quantum in- 
formation theory the physical conditions underlying the 
Markovian approximation can be violated since the cou- 
pling of the reduced system to the environment may be 
of the same order of magnitude as the environment re- 
laxation rate. Consequently, the time evolution of the re- 
duced system is no longer Markovian and the equations 
of motion either become temporally non-local Nakajima- 
Zwanzig equations [13, [13 oi' explicitly time dependent 
time-convolutionless equations [Tjl , depending on the ap- 
proach chosen 'l'|. A feature that these equations share 
is that they are often difficult to derive for any given 
model, and also to solve. Moreover, they may lead to 
non-physical behavior such as the violation of positivity 
of the dynamical map [3, [la, I16, 17, 18, 19] which is 
a consequence of the phenomenological nature of most 
of the non-Markovian approaches. Since there are cur- 
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FIG. 1: Model for the non-Markovian decay of an entangled 
two-atom system. Two (entangled) atoms 1 and 2 are indi- 
vidually coupled to two additional atoms 3 and 4, which are 
in turn coupled to a thermal reservoir with decay rates 7. 



rently no established criteria that can be used to identify 
non-Markovian master equations that preserve complete 
positivity, the only way to ensure this property is to find 
a model that itself generates completely positive maps 
and to use it with no approximations in order to derive 
the final equations of motion. This has recently been 
achieved by Salo, Barnett and Stenholm [23 by consid- 
ering the problem of non-Markovian thermal damping 
of a two-level atom. The model consisted of coupling 
the principal atom to an additional two-level atom which 
in turn is coupled to a thermal environment. Utilizing 
the exact Nakajima-Zwanzig elimination procedure en- 
ables them to find a master equation with memory that 
is guaranteed to generate completely positive evolution. 
In this paper we discuss the dynamics of the entan- 
glement of a bipartite system in a non-Markovian reser- 
voir. The system is composed of two qubits which are 
represented by two- level atoms. In order to derive an ex- 
act non-Markovian master equation we extend the model 
of Salo et al. 20] to an entangled two-atom system to 
guarantee its completely positive evolution. As in [20l] . 
the two atoms are individually coupled to two additional 
two-level atoms which in turn interact individually with 
a thermal environment, as shown in Fig.^ This ensures 
that the reservoirs are not able to entangle the two-atom 



system and, consequently, it is possible to derive non- 
Markovian master equations which are of a similar form 
to the equations in j2.(|. However, although the reservoir 
is not able to entangle the two-atom systems the dynam- 
ics of the entanglement of an initially entangled two-atom 
system shows a behavior which clearly demonstrates the 
memory effects of the non-Markovian reservoir. In par- 
ticular, we show that the entanglement, quantified by 
the concurrence [22, 123 , shows an oscillating behavior in 
time that depends on the degree to which the reservoir 
is non-Markovian. Thus, the reservoir remembers that 
the two-atom system was initially entangled, and the en- 
tanglement, which may vanish during the time evolution, 
can be retrieved by the system, leading to a "collapse and 
revival like" time evolution of the entanglement. Such 
an behavior may be of importance in quantum compu- 
tational systems that are exposed to non-Markovian de- 
coherence effects because even the concurrence of an ini- 
tially completely entangled state can vanish after a finite 
period of time [23, |2J|. Note, however, that the con- 
currence will partially revive due to the memory of the 
reservoir that feeds back into the system. Note further, 
that this model allows decoherence effects to be modi- 
fied by adjusting the parameters between the interacting 
atoms. This is to some extent similar to the modification 
of decoherence |2^ |2^ in engineered reservoirs [23, llg • 
The rest of this paper is organized as follows: In Sec. 
^we introduce the parameterized thermalization model 
for the entangled qubit system. In Sec. lIIII we derive the 
dynamical equations of the coupled two-atom system and 
solve them numerically, assuming that the system is ini- 
tially in a maximally entangled quantum state of the Bell 
form. We discuss the dynamical evolution of the entan- 
glement which is quantified by the concurrence, in Sec. 
IIVI We also discuss the effects of a finite-temperature 
non-Markovian reservoir on the dynamical evolution of 
the concurrence. We conclude in Sec. with some dis- 
cussion. 



II. PARAMETERIZED THERMALIZATION 
MODEL FOR TWO ENTANGLED QUBITS 

In this section we derive a simple toy model consisting 
of two entangled qubits that are individually coupled to 
non-Markovian thermal reservoirs. This model is based 
on a recently proposed system that consists of a single 
two-level atom coupled to a specific non-Markovian ther- 
mal bath, such that it is possible to derive an exact mas- 
ter equation that preserves complete positivity |20j . The 
system is composed of a two-level atom 1 that is coupled 
via an exchange interaction to another atom 2, which, in 
turn, is damped by a thermal reservoir at a specific tem- 
perature that can be characterized with the bosonic exci- 
tation number n. We extend this model to an entangled 
two-atom system 1 and 2 where both atoms are assumed 
to be two-level atoms which are individually coupled to 
two additional and possibly different two-level atoms 3 



and 4 via an exchange interaction. The atoms 3 and 4 
are both damped by a thermal reservoir at a temperature 
that is characterized with the bosonic excitation number 
n. Consequently, the atoms 3 and 4 play the role of a 
"near-environment" to the atoms 1 and 2, respectively. 
Their dynamical influence produces memory effects in 
the subsystem of atom 1 and 2. The individual coupling 
of the atoms to the reservoir ensures that the dynamics 
of the two atoms can be treated independently. We note, 
that in the single-atom case at zero temperature with no 
initial environment photons present, this system corre- 
sponds to the damped Jaynes-Cummings model [29J, as 
discussed in [20l| . 

Suppose the atoms 1 and 2 are initially in an arbitrary 
state. Further, assume that the atoms 3 and 4 are ini- 
tially uncorrelated with the atoms 1 and 2 and that they 
are both in a thermal state 
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The initial state of the total system is thus described as 

p{t == 0) = pi234(t = 0) = pi2(t = 0) P3 P4. (2.2) 

Here, pi2 is an arbitrary density operator of the subsys- 
tem 1 — 2 that allows us to describe initially entangled 
states of atoms 1 and 2. The atomic states represent 
the computational qubit basis |0) and |1) for the atomic 
ground and excited state, respectively. 

The dynamics of the composite four-atom system is 
governed by the following equation (h = 1) 
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dt 
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H2 + Hi + H^!t\pit)] + Ap(i), (2.5) 
c^„|l)„„(l|, n= 1,2,3,4 (2.6) 

"1 (^^i^cTs" + ^i^^s) ' (2-7) 

(2.8) 



= a2 (o'^Cr4 + 0-2 CT 



Csp ^ 7 (n + 1) (2(73^ pcr+ -cr+CT^p-pcr+CT^) 

+jn{2a^pa^ - cr^a+p - p(T^a+), (2.9) 

C4P = "f{n+l){2al pa'l - a'lal p- patcrl) 

+^n{2a4 pa^ - a^a^p - pa^at)- (2.10) 

Here, w„, n = 1, 2, 3, 4 is the energy of the excited state 
of atom n, ai and a2 are the coupling strengths between 
atoms 1 — 3 and 2 — 4, respectively, 7 is the thermaliza- 
tion rate for atoms 3 and 4, and n is the thermal bosonic 



excitation number. In order to simplify matters we have 
assumed that the atoms 3 and 4 are coupled to a sin- 
gle reservoir, although in principle they could be coupled 
to separate reservoirs with different temperatures. Note, 
that the dynamical equation (|2.3I) is unable to generate 
entanglement between the atoms 1 and 2 since the equa- 
tions of motion decouple between the atomic subsystems 
1 — 3 and 2 — 4. However, if entanglement is initially 
present in the atomic subsystem 1 — 2, the memory effects 
which emerge from the coupling to the atomic subsys- 
tems 3 and 4 will be able to modify the decoherence rate 
of the entanglement. Consequently, the non-Mar kovian 
bath dynamically influences the decoherence rate of the 
entanglement of subsystem 1 — 2 because of the feedback 
of the atomic subsystems 3 and 4 on the system 1 — 2 
under consideration. 

The underlying idea of this non-Markovian bath model 
is to assume the atoms 3 and 4 serve as "mini-reservoirs" 
which arc able to dynamically modify the equations of 
motion of the relevant atoms 1 and 2. The coupling 
between the mini-reservoirs can be steered via the in- 
teraction parameters ai and a2 as well as the detuning 
parameters Afc — Wk — uJk+2 for k = 1,2. These mini- 
reservoirs will be unobserved and because of their cou- 
plings to a thermal Markovian reservoir they are assumed 
to remain in a certain reference state, namely a thermal 
reservoir state of a temperature that is imposed by the 
thermal environment. Their dynamical influence which 
feeds back into the subsystem 1 — 2, however, introduces 
non-Markovian effects in the atomic subsystem 1 — 2. 
This remains true even in the case when they remain un- 
observed and are assumed to be in the thermal reference 
state, as it has been demonstrated in |20l |. 

The Nakajima-Zwanzig equation of motion \VA ll^ 
projects the combined quantum system into the rele- 
vant parts, described by atom 1 and 2, and the irrelevant 
parts, represented by atoms 3 and 4, utilizing projector 
operators Vk and Qfc, where fc = 1, 2 indicate the atomic 
subsystems 1 and 2. When we insert the thermal state as 
the reference state for the atoms 3 and 4, the projectors 
are given as 



Vipiz 


= Tr3[pi3](g)p3, 


(2.11) 


'P2P2i 


= Tr4[p24]«)p4, 


(2.12) 


QlPl3 


= Pl3 " Tl-slPls] <»P3, 


(2.13) 


Q2P24 


= (024 - Tr4[p24] <8)P4, 


(2.14) 



where pia and p24 define arbitrary operators (not neces- 
sarily density operators) of the atoms 1 — 3 and 2 — 4, 
respectively. Note, that in general, the projectors can 
operate on arbitrary operators of the combined systems 
1 — 3 and 2 — 4. The Nakajima-Zwanzig equation of mo- 
tion has been derived for the single atom case in |20| 
under the assumption that, initially, the atom is not en- 
tangled with the "mini-reservoir atom" at i = and the 
mini-reservoir atom is in the thermal state p. In the two- 
atom case, the situation which is discussed in this paper, 
the Nakajima-Zwanzig equations of motion turn out to 



be of the same form. 
d 
di 



VkP^^Ht) = VkC^'^VkP^'Ht) 



+ / Pfc/:^'^a('Hi-r)Qfc/:(^¥fep('=)(r)dr. 
Jo 

(2.15) 

This follows from the fact that the dynamical equation 
of motion H2.3|l is not able to entangle the atoms 1 and 2. 
Consequently, we can rewrite every density operator of an 
entangled state between atoms 1 and 2 as a linear combi- 
nation of tensor products of operators on systems 1 and 2 
for which the single atom Nakajima-Zwanzig formalism 
of _20] applies. Here, g('=)(i - r) = exp[(i - T)Qfc£('=)] 
where it is assumed that C^'^' is time-independent and 
p'^^^t)^ fc = 1,2 represents the operators pi^ and p24, 
respectively. 

The Nakajima-Zwanzig equation of motion has been 
derived in |20l| to which we refer for further details. What 
is important in these master equations is that the oper- 
ator expressions describing the memory terms are not 
of the Lindblad form. In addition, the memory kernel 
contains two different memory functions one for the di- 
agonal part of the density operators and another one for 
the coherences. Hence the total convolution kernel can 
not be of the form K{t ~ t)C where £ is a Lindblad op- 
erator. Nevertheless, under the assumption (|2.2I) of the 
initial state of the total combined system the Nakajima- 
Zwanzig equation of motions is an exact equation that 
can be derived without any approximations and its so- 
lutions remain legitimate density operators at all times. 
In other words the dynamical equation of the subsystem 
1 — 2 is described by a completely positive map. In the 
next section we solve the master equation (|2.3ll by pro- 
jecting the total density operator onto the subspace in 
which the atoms 3 and 4 remain in the thermal reference 
state given in Eq. I|2.1|) during time evolution. This en- 
sures that the density operator puit) of the subsystem 
1 — 2 is a legitimate density operator at all times. 



III. DYNAMICAL EQUATIONS OF THE 
COMPOSITE TWO-ATOM SYSTEM 

In this section we numerically solve the dynamical 
equations of the density operator pi2(t) of an initially 
entangled two-atom system in the non-Markovian reser- 
voir which is represented by the model of the preceding 
section. Initially, we suppose the two atom system to be 
in a completely entangled Bell state. 
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The corresponding density matrix is accordingly given as 
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for all k e 1,2. We represent the density matrix of the 
total system pi234(i) using the operator basis derived in 
|2(ll |. In this basis we can express 
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where the operator basis satisfies, 
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Clearly, in the representation (|3.4I) above, the dynamical 
equations H2.3|) decouple in systems (1) and (2) and we 
describe the effect of the generator of motion O''^ on the 
systems A: = 1, 2, as 
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Here, we have introduced the effective decoherence rate 
7cff = (2n+ 1)7 as well as the detuning A^ — ujk — ujk+2- 
The matrix representation of the generator £"'^ in the 

operator basis < Xm \ is thus 
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The dynamical equation of the density operator H3.4|) follows from H2.3|l : 
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and it decouples for fc = 1 and fc = 2, so expressed as a vector, the coefficients 
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\t),---,clt\t) 
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are determined by the differential equation 



dt 



:('=)(t) = £('=) cW(i), 
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where C^'''> is given in H3.7(l . The relevant part of 
the equation of motion of system (fc) is represented in 

the subspace spanned by X^ ,X[ t-^5 ' ^^'^ -^7 : 
i.e. by the dynamical evolution of the coefficients 
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{t),c\ ' (i), C5 ' (i), and C7 ' (t). Consequently, the pro- 
jection operators V'^^^ and QP^^ are defined as 
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where (fc) indicate that the projection operators oper- 
ate exclusively on the subsystems (fc). The equation of 
motion of the relevant density operator 

Pl234(t) = V'^^^V'^^'' pi234it) = Tr34[pi234(t)] (» P3 '^P4, 

(3.13) 
therefore describes the subspace in which the mini- 
reservoirs represented by atoms 3 and 4 remain in the 
thermal reference state p at all times. Note, however, 
that the dynamical evolution of the coefficients repre- 
senting the relevant part of the dynamics is influenced 
by the dynamical evolution of coefficients that represent 



the irrelevant part of the dynamics. This ensures that 
memory effects of the reservoir are imposed on the dy- 
namics of the relevant subsystem, i.e. the reservoir feeds 
back into the subsystem. 

The initial state of the total system pi234(i ~ 0) is 
given as 



Pl234{t = 0) = pi2(t = 0) «) P3 (g) P4, 



(3.14) 



where, pi2{t = 0) is defined in (|3.2|l . Since the differen- 
tial equation of motion is linear, we can solve it for each 
tensor product term in (|3.2|) independently, and add up 
the solutions. For the first tensor product contribution 
|l)i(l| (X) 11)2(11 of pi2{t = 0) we arrive at the vector of 
coefficients 
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(3.15) 
For the second contribution |0)i(0| |0)2(0|, we obtain 
the following initial condition 
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lowing coefficients 
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)2(1| leads to the fol- 



f(i)(t = 0) = [0,0,0,0,0,0,0,1,0]' 
while for the last contribution |l)i(0| 

r(l) 



= f(2)(t = 0), 

(3.17) 
11)2(01 we obtain 

(2) 



g^^^(i = 0) = [0,0, 0,0, 0,1, 0,0,0]^ =g(2)(t^0). 

(3.18) 
The time-evolution of the relevant density operator in the 
appropriate subspace projected to P3 (g) P4 is described 
by the dynamical evolution of the projected vectors 
of coefficients pWcW(i),p('=)dW(i),pWfW(i), and 
7:'(fc)g(*:)(i)^ where fc = 1, 2. We calculate the time evolu- 
tion of the vectors of coefficients c^'^\t),d^''^{t),f^''\t), 
and g'-'^^(i) numerically from the equation of motion 
(|3.10l) . with the initial conditions given in H3.15|l - 13.18ll . 
The time evolution of the relevant density operator (|3.13() 
is then given by 
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FIG. 2: Time evolution of the density matrix elements for the 
parameters Ai = A2 = 2, ai = a2 = 2, 7 = 0.5 and n — 0. 
The frequency wi is assumed to be 10. The structure of the 
density matrix remains of the same form given in Eq. (14.61 1 
for all parameters used in this paper. 



This describes the time evolution of the density opera- 
tor (|3.4() projected onto the subspace where the atoms 3 
and 4 remain in the thermal state H2.1|l at all times. The 
evolution of the reduced systems operator pi2 (t) of the 
subsystem 1 — 2 follows from Eq. (|3.19() by tracing out the 
near reservoir systems 3 and 4. The dynamical evolution 
of the reduced density operator pi2{t) = Tr34[pi234(i)] in 
matrix representation of the computational basis is dis- 
played in Fig.|5]for the initial condition given in l|3.2() and 
for certain interaction parameters with the near environ- 
ments. Of course, pi2it) describes a real density matrix 
that satisfies the positivity condition since the map from 
which it has been derived is completely positive. 



IV. DYNAMICS OF ENTANGLEMENT 
QUANTIFIED BY CONCURRENCE 

In this section we consider the dynamics of the en- 
tanglement of the bipartite system 1 — 2 in the non- 
Markovian reservoir assuming the system density opera- 
tor pi2{t = 0) to be in the state H3.2|l initially. We take 
Wootters concurrence [2ll [23 to quantify the amount of 
entanglement present in the subsystem pi2(t). The con- 



currence of a bipartite, arbitrary density operator p of 
two qubits is defined as 



C{p) = max {0, Ai - A2 - A3 - A4} , 



(4.1) 



where the Ai are the square roots of the eigenvalues of 
p/5 in descending order. Here p is the result of apply- 
ing the spin-flip operation to the complex conjugation or 
transpose of p: 



P^ {(7ylg) ay)p*{(Jy ® (7y), 



(4.2) 



where the complex conjugation is taken in the standard 
computational basis. The concurrence is related to the 
entanglement of formation Ef{p) by the following func- 
tion i21j,|2^ 



where, 
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and 



h{x) = — xlog2 X — {\ ~ x) log2(l — x). 
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The structure of the density matrix of the two-atom 
system remains in the following form during the time- 
evolution 
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c{t) 
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where a(i) + 6(t) + c(f) + rf(f) = 1. This property is a con- 
sequence of the special initial condition and the fact that 
the environment, regardless of whether it is Markovian 
or not, can not create coherence that was not present ini- 
tially. The time evolution of the density matrix is shown 
in Fig. |2] for a certain set of parameters. This particular 
form of the density matrix (|4.(j|| allows us to analytically 
express the concurrence as 



C[pi2(t)] = 2max{0, |/(t)| - ^h{t)c{t)}. (4.7) 

In Figs. I3I()I we plot the evolution of the concurrence 
for different parameters of the non-Markovian bath. As 



a result of the memory effects the concurrence displays 
an oscillatory time dependence, especially when the cou- 
pling parameter a is strong. In the case of zero detun- 
ing, A = 0, the entanglement can vanish after a finite 
time even when the initial state of the two-atom system 
is a maximally entangled Bell state, as shown in Fig. 
01 However, due to the memory inherent in the non- 
Markovian bath, the entanglement is partially restored. 
The two-atom system remembers its initial degree of en- 
tanglement to some extent, which leads to a revival of 
the entanglement. The collapse and revival of the entan- 
glement displays a damped oscillatory behavior in time, 
whose frequency depends on the coupling parameter a. 
Thus the coupling parameter a characterizes the memory 
ability of the reservoir to some extent. In case when the 
detuning A differs from zero, the entanglement does not 
completely vanish after a finite time, as displayed in Fig. 
13 The detuning suppresses the decoherence and thus 
extends the decoherence time. This is because the mini- 
reservoir atoms 'shield' atoms 1 and 2 from the reser- 
voir. When the damping rate 7 is of the same order 
or greater than the coupling parameter a, the decay of 
the entanglement is damping oriented and oscillations do 
not appear, as seen in Fig. |S1 However, the effect of the 
non-Markovian bath is still visible, as the concurrence 
displays a non-exponential decay in time, and the deco- 
herence time is extended. The last plot Fig. |H| displays 
the concurrence in the Markovian limit, which is achieved 
when a/7 — > but a^/7 remains finite 20] . In this case 
the Markovian time scale is given as 20J 

1 „.2 

(4.8) 



(2n+l)2 7 ' 

leading to an effective decoherence rate of F = 1/3 for 
the parameters used in Fig. El 

An important observation with respect to entangle- 
ment is that for an open quantum system, entanglement 
can die off at finite times. It has recently been shown that 
this can happen not only when one starts with mixtures 
of entangled quantum states, but also in cases of initially 
pure and maximally entangled states, if the entangled 
system is coupled to a thermal reservoir at finite temper- 
ature [23, |2J, I33 • In other words non-local disentangle- 
ment times are in general shorter than local decoherence 
times. In practical quantum computation schemes which 
operate at finite temperature this is of particular impor- 
tance since some quantum computational algorithms rely 
on entanglement. We therefore ask next, to what ex- 
tent a non-Markovian reservoir affects the decoherence 
of the entangled system in the case of an environment 
with non-zero temperature. This is illustrated in Figs. 
1711 II for an environment with a mean photon number of 
n = 0.2, i.e. a hot reservoir. In order to visualize the 
periods over which the concurrence becomes zero more 
clearly, we do not display the actual concurrence given 
in gTj) but the quantity 2{\f{t)\ - [b{t)c{t)]^/'^}, which 
is less than zero when the concurrence vanishes. In con- 
trast to the zero temperature case we can see that the 




FIG. 3: Evolution of the concurrence of an initially maximally 
entangled Bell state which decays into the non-Markovian 
reservoir. The dimensionless parameters of the reservoir are 
Ai = A2 = 2, ai = 02 = 2, 7 = 0.5 and n = 0. The 
frequency ui is assumed to be 10. The concurrence displays 
and oscillating behavior due to the memory effects of the non- 
Markovian reservoir. Entanglement which is initially lost due 
to decoherence is partially restored by the memory effects. 
Because of the large detuning, entanglement does not com- 
pletely vanish in a finite time. 
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FIG. 4: Evolution of the concurrence of an initially maximally 
entangled Bell state which decays into the non-Markovian 
reservoir. The dimensionless parameters of the reservoir are 
Ai = A2 = 0, Qi = 0:2 = 2, 7 = 0.5 and n = 0. The fre- 
quency ui is assumed to be 10. As in Fig. |21the concurrence 
displays and oscillating behavior due to the memory effects 
of the non-Markovian reservoir. However, the entanglement 
completely vanishes at finite times before it is partially re- 
covered in contrast to Fig. |21 This is a consequence of the 
vanishing detuning. 



concurrence vanishes suddenly after a finite time for all 
plots, a behavior which has been predicted for Marko- 
vian reservoirs in |23, 24] . In addition, in cases of strong 
coupling to the near environments, when a > 7, the 
concurrence can vanish for a period of time, before be- 
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FIG. 5: Evolution of the concurrence of an initially maximally 
entangled Bell state which decays into the non-Markovian 
reservoir. The dimensionless parameters of the reservoir are 
Ai = A2 = 0, ai = Q2 = 0.5, 7 = 0.5 and n = 0. The 
frequency uji is assumed to be 10. The concurrence does not 
display an oscillating behavior any more since the interaction 
with the near environment is comparable weak and damping 
dominates. However, the decay of the concurrence is non- 
exponential which is a manifestation of the non-Markovian 
reservoir. 
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FIG. 6: Evolution of the concurrence of an initially maximally 
entangled Bell state which decays into a Markovian reservoir 
which is achieved when a/7 — > 0. The dimensionless param- 
eters of the reservoir are Ai — A2 = 0, ai = 0:2 = 3, 7 = 27 
and n = giving an effective Markovian decay rate of F = 1/3 
(see text). The frequency wi is assumed to be 10. 
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FIG. 7: Evolution of the concurrence of an initially maximally 
entangled Bell state which decays into the non-Markovian 
reservoir of finite temperature with mean photon number 
n = 0.2. Here, we do not plot the actual concurrence 
but C[/9i2(t)] = 2{|/| - [b(t)c(f)]i/2| instead of C[pi2{t)] = 
2max{|/| — [b{t)c{t)Y'^} > to see how negative it be- 
comes. The dimensionless parameters of the reservoir are 
Ai = A2 = 0, ai = a2 = 5, 7 = 1/3 and n = 0.2. The 
frequency loi is assumed to be 10. As in Figs. |3 and 0] the 
concurrence displays and oscillating behavior due to the mem- 
ory effects of the non-Markovian reservoir. In contrast to the 
zero temperature case, however, the concurrence vanishes at 
a finite time completely and cannot be retrieved. In addition, 
during the time period where the concurrence can be partially 
recovered there are extended periods of time where it vanishes 
completely. 




FIG. 8: Same as Fig. |7|but with parameters Ai — A2 = 2, 
ai = 02 = 2, 7 = 0.5 and n = 0.2. The frequency ui is 
assumed to be 10. 



ing partially recovered due to the memory effect of the 
non-Markovian reservoir. As the interaction parameter 
a with the near environment becomes stronger, these pe- 
riods of time when concurrence completely vanishes be- 
come shorter, and the rate at which the concurrence par- 
tially revives increases as shown in Figs. I7I9I The time 
until there is no more revival of the concurrence becomes 
longer when there is some detuning between the system 
under consideration and the near environment as can be 
seen by comparison of Figs. |H1 and |^ When 7 > a and 
the damping dominates, the effect of the non-Markovian 



reservoir becomes weaker and the evolution of the concur- 
rence does not display an oscillating behavior any more. 
However the decay of the concurrence is non-exponential 
and the period of time before it vanishes completely can 
be extended in comparison with a Markovian reservoir, 
as can be seen by comparing Fig. ^| with Fig. lllL which 
shows the Markovian limit. For the parameters given 
in Fig. E] the effective Markovian decay rate F for the 
Markovian limit can be calculated from Eq. (|4.8|) to be 
F w 1/6. Although this effective decay rate is consid- 





FIG. 9: Same as Fig. [7| but with parameters Ai = A2 = 0, 
ai = a2 = 2, 7 = 0.5 and n = 0.2. The frequency wi is 
assumed to be 10. 
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FIG. 10: Same as Fig. |7|but with parameters Ai = A2 = 0, 
ai = a2 = 0.5, 7 = 1 and n = 0.2. The frequency ui is 
assumed to be 10. This figure displays the strong damping 
case in which the efi'ect of the near environment is comparably 
weak. However non-Markovian effects are clearly visible since 
the decay of the concurrence is non-exponential. As a result 
the period of time during which the concurrence is maintained 
is extended compared to the Markovian decay. 



erably smaller than the decay rate 7 = 1 used in the 
non-Markovian reservoir of Fig. IIUI the period of time 
before the concurrence suddenly vanishes is shorter than 
for the non-Markovian reservoir. Thus a non-Markovian 
reservoir helps to maintain concurrence over a larger pe- 
riod of time in comparison to a Markovian reservoir. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have discussed the effects of a non- 
Markovian bath on the decoherence of an entangled two- 
atom system where the atoms are assumed to be two-level 
systems representing qubits. The non-Markovian bath is 
described by a model that leads to an exact Nakajima- 
Zwanzig type master equation that can be solved with- 
out any approximations, as discussed by Salo et al. in 



FIG. 11: Evolution of the concurrence in the Markovian limit 
of the temperature reservoir. The parameters are Ai = A2 = 
0, Qi = a2 = 3, 7 = 27 and n = 0.2 leading to an effective 
Markovian decay rate of F ~ 1/6. The frequency lji is as- 
sumed to be 10. In comparison to Fig. llOl the concurrence is 
maintained during a shorter period of time. 



|2nl |. The advantage of this master equation is that it 
leads to a completely positive time evolution and there- 
fore non-physical quantum states do not arise, a situa- 
tion which can appear in parameterized phenomenologi- 
cal non-Markovian master equations [l^ . 

We discuss the time-evolution of the entanglement of 
an initially maximally entangled Bell state in this non- 
Markovian reservoir. The entanglement is quantified by 
the concurrence |2ll |22| which can be derived in an an- 
alytical form since the density matrix of the two-atom 
system retains a certain structure during time evolution. 
The dynamics of the concurrence displays some inter- 
esting features which are manifestations of the memory 
effects inherent in a non-Markovian reservoir. In par- 
ticular, the entanglement can vanish after a finite time 
but, depending on the parameters of the reservoir, it can 
be partially restored. The revival of the concurrence is 
an effect that is specific to the non-Markovian reservoir 
since the reservoir remembers that the system was ini- 
tially completely entangled, and some of this entangle- 
ment can be restored in the principle system even though 
it becomes zero during the time evolution. 

In the case of finite temperature baths it is well known 
that the concurrence of an entangled system can expe- 
rience "sudden death" when the system decays into a 
Markovian reservoir [23, Hj, ISfl . This sudden death of 
the concurrence can not be prevented in a thermal non- 
Markovian reservoir either. However, due to the memory 
effect in the non-Markovian reservoir the concurrence can 
revive even when it completely vanished after a finite pe- 
riod of time. Depending on the parameters used, the con- 
currence displays a collapse and revival structure in time, 
which finally ends in sudden death. The time of the sud- 
den death can be extended by the non-Markovian reser- 
voir in comparison with the Markovian reservoir. This is 
an effect of the non-exponential decay of the concurrence 
in a non-Markovian reservoir compared with its exponen- 
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tial decay in a Markovian reservoir. 
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